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Quantum Monte Carlo (QMC) and Maximum Entropy (ME) techniques are used to study the 
spectral function A(p, uj) of the one band Hubbard model in strong coupling including a next- 
nearest- neighbor electronic hopping with amplitude if /t = —0.35. These values of parameters 
are chosen to improve the comparison of the Hubbard model with angle-resolved photoemission 
(ARPES) data for Sr2Cu02Cl2. A narrow quasiparticle (q.p.) band is observed in the QMC 
analysis at the temperature of the simulation T = t/3, both at and away from half-filling. Such 
a narrow band produces a large accumulation of weight in the density of states at the top of the 
valence band. As the electronic density (n) decreases further away from half-filling, the chemical 
potential travels through this energy window with a large number of states, and by (n) ~ 0.70 
it has crossed it entirely. The region near momentum (0, 71") and (it, 0) in the spectral function is 
more sensitive to doping than momenta along the diagonal from (0, 0) to (n, tt). The evolution with 
hole density of the quasiparticle dispersion contains some of the features observed in recent ARPES 
data in the underdoped regime. For sufficiently large hole densities the "flat" bands at (it, 0) cross 
the Fermi energy, a prediction that could be tested with ARPES techniques applied to overdoped 
cuprates. The population of the q.p. band introduces a hidden density in the system which produces 
interesting consequences when the quasiparticles are assumed to interact through antiferromagnetic 
fluctuations and studied with the BCS gap equation formalism. In particular, a region of extended 
s-wave is found to compete with d-wave in the overdoped regime, i.e. when the chemical potential 
has almost entirely crossed the q.p. band as (n) is reduced. The present study also shows that 
previous "real-space" pairing theories for the cuprates, such as the Antiferromagnetic Van Hove 
scenario, originally constructed based on information gathered at half-filling, do not change their 
predictions if hole dispersions resembling non-interacting electrons with renormalized parameters 
are used. 

PACS numbers: 74.20.-z, 74.20.Mn, 74.25.Dw 



I. INTRODUCTION 

Recent studies of E^S^Cai-xDyxC^Os+a using 
angle-resolved photoemission techniques have provided 
the evolution of the quasiparticle band as the hole density 
changes in the underdoped regime of the high-T c cuprates 
|0,Hj3[| . These studies complement previous ARPES anal- 
ysis of the hole dispersion in the antiferromagnetic insula- 
tor Sr2Cu02Cl2 IHJ^]. The overall results emerging from 
these experiments can be summarized as follows: (i) the 
bandwidth of the quasiparticle band is a fraction of eV, 
i.e. narrower than generally predicted by band structure 
calculations. This result suggests that strong correla- 
tions are important in the generation of the quasiparticle 
dispersion; (ii) the two dimensional (2D) t — J model 
explains accurately the bandwidth for a hole injected in 
the insulator, as well as the details of the dispersion along 
the main diagonal in momentum space from p = (0, 0) 
to (tt,tt) (in the 2D square lattice notation) Q; (iii) how- 



ever, the behavior along the direction from (0, tt) to (tt, 0) 
is not properly described by the t — J model which pre- 
dicts a near degeneracy between momenta along this line 
||, in contradiction with experiments. The ARPES data 
for the antiferromagnetic insulator show that the quasi- 
particle signal at (tt, 0) is very weak, and about 0.3 eV 
deeper in energy than (7r/2,7r/2) j|; (iv) as the den- 
sity of holes grows the region near (tt, 0) moves towards 
the chemical potential which is approximately reached 
at the optimal concentration while the main diago- 
nal (0,0) — (ir,ir) is less affected; (v) at optimal doping 
remarkably "flat" bands at (tt, 0) near the Fermi energy 
have been observed 0. 

The results of Marshall et al. jl| have been interpreted 
as providing evidence for "hole pockets" in the under- 
doped regime caused by short-distance antiferromagnetic 
correlations, although other explanations such as pre- 
formed pairs are also possible ^||. As explained be- 
fore, the standard t — J model is not enough to provide 
all the details of the ARPES results at half-filling, and 
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thus presumably it cannot explain the data away from 
half-filling either. However, the t — J model is just one 
possible Hamiltonian to describe the behavior of holes 
in an antiferromagnctic background. While the model 
is attractive for its simplicity, there is no symmetry or 
renormalizability argument signaling it as unique for the 
description of the cuprates. For this reason Nazarenko 
et al. JlO| recently introduced extra terms in the t — J 
Hamiltonian to improve the agreement with experiments. 
Adding an electronic hopping along the diagonals of the 
elementary plaquettes (with amplitude t' jt — —0.35) the 
results for the insulator were improved since the posi- 
tion of the quasiparticle at (tt, 0) proved to be very sen- 
sitive to the strength of extra hole hopping terms in the 
model. Actually, the q.p. peak at this momentum moves 
towards larger binding energies as a negative t' /t grows 
in amplitude ju|. This idea has been used by several 
other groups which, in addition to hoppings along the 
plaquette diagonals, have included hoppings at distance 
of two lattice spacings with amplitude t" to obtain an 



even better agreement with experiments |11 



The purpose of this paper is to report on results for 
the density dependence of the quasiparticle dispersion 
corresponding to the one band Hubbard model with a 
next-nearest-neighbor hopping t' working in the strong 
Coulomb coupling regime. It is expected that the U —t—t 1 
and t — t' — J models produce qualitatively similar physics 
at large U ft. As numerical technique we use the Quan- 
tum Monte Carlo method supplemented by Maximum 
Entropy analysis. The success of the t — t' — J model at 
half-filling p)||il"| and the recent availability of ARPES 
data in the underdoped regime |U|§|| prompted us to 
carry out this study. Note that although the numerical 
methods used here are powerful, their accuracy is lim- 
ited and, thus, our results are mostly qualitative rather 
than quantitative. Nevertheless, from the analysis of the 
present QMC data and also comparing our results against 
those produced previously for the t — J and Hubbard 
models here we arrive to the conclusion that at least 
some features of the experimental ARPES evolution of 
the hole dispersion can be explained using one band elec- 
tronic models. Experimental predictions are made to test 
the calculations. The presence of a narrow quasiparticle 
band in the spectrum, and its crossing by the chemical 
potential as the electronic density is reduced, are key fea- 
tures for the discussions below. 

The present paper is not only devoted to the analysis of 
QMC results and its comparison with ARPES data, but 
it also addresses the influence of hole doping in recently 
proposed theories for the cuprates jl2| . In these theories 
the hole dispersion is calculated at half-filling and as- 
sumed to change only slightly as the hole density grows. 
Hole attraction is assumed to be dominated by the mini- 
mization of antiferromagnetic (AF) broken bonds, which 
implies the presence of an effective nearest-neighbor (NN) 
density-density attraction proportional to the exchange 
J. Superconductivity in the d x 2_ y 2 channel is natural in 
this scenario |L2| due to the strong AF correlations P-fbl . 



The interaction of holes is better visualized in real space, 
i.e. with pairing occurring when dressed holes share a 
spin polaronic cloud, as in the spin-bag mechanism Jl4| . 
This real-space picture (see also Ref. |15|]) holds even for 
a small AF correlation length, £,af, and in this scenario 
there is no need to tune parameters to work very close to 
an AF instability as in other approaches. 

To obtain quantitative information from these intuitive 
ideas, holes moving with a dispersion calculated using 
one hole in an AF background, €af(jp), and interacting 
through the NN attractive potential mentioned before, 
have been previously analyzed fHq] . Within a rigid band 
filling of €af(p) and using a BCS formalism, d x 2_ y 2 su- 
perconductivity dominates with T c ~ 100K caused by a 
large density of states (DOS) that appears in the hole 
dispersion. The idea has many similarities with previous 
scenarios that used van Hove (vH) singularities in the 
band structure to increase T c [na, although d-wave su- 
perconductivity is not natural in this context unless AF 
correlations are included. However, the rigid band fill- 
ing is an approximation whose accuracy remains to be 
tested. In particular, the following questions naturally 
arise: (i) does the q.p. peak in the DOS found at (n) = 1 
survive a finite hole density?; (ii) to what extent do the 
changes in the q.p. dispersion with doping affect previous 
calculations in this framework?; (iii) are the "shadow" re- 
gions generated by AF corr elations @,|l|,|l|,|l9| impor- 
tant for real-space pairing approaches? With the help 
of the present QMC results, as well as previous simula- 
tions, in this paper all these issues are discussed. The 
overall conclusion is that as long as the hole density is 
such that the £af is at least of a couple of lattice spac- 
ings, the predictions of previous scenarios |T^] and other 
similar theories remain qualitatively the same, in spite of 
substantial changes occurring in the hole dispersion with 
doping. 

The organization of the paper is as follows: in Sec. II 
the model and details of the numerical method are dis- 
cussed. Sec. Ill contains the QMC-ME results both at 
half- filling and with a finite hole density. In Sec. IV the 
results are discussed and compared with ARPES experi- 
ments. Implications for real-space theories of high-T c are 
extensively studied in Sec.V. Results of previous publi- 
cations are also used to construct a simple picture for 
the behavior of electrons in the Hubbard model. Sec. VI 
contains a summary of the paper and its experimental 
predictions. 



II. MODEL AND NUMERICAL TECHNIQUE 

The one band Hubbard Hamiltonian with next- 
nearest-neighbor hopping is given by 

h = - 1 (4,o c s,° + h - c -) - *' Yl ( c L c s<° + /i - c -) 

<ij>,^ «ij»,<7 
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-C/5> iT -l/2)(n u -l/2) + M £<, 



(1) 



where cj , 



creates an electron at site i with spin pro- 



jection <t, riin is the number operator, the sum (ij) runs 
over pairs of nearest-neighbor lattice sites, ((ij)) runs over 
pairs of lattice sites along the plaquette diagonals, U is 
the on-site Coulombic repulsion, t the nearest-neighbor 
hopping amplitude, t' is the plaquette diagonal hopping 
amplitude, and /x the chemical potential. Throughout 
this study we will set t = 1, t'/t = —0.35, and use peri- 
odic boundary conditions. 



Using standard QMC methods |20|, we have obtained 
the imaginary time Green's functions at finite temper- 
ature. The method of Maximum Entropy pl| was used 
to analytically continue the imaginary time Green's func- 
tions to obtain the spectral weight function A(p, uj). Pre- 
vious studies using this method have concentrated on, 
e.g., the one-dimensional ]22| ] and two-dimensional Hub- 
bard model at and away from half-filling [T^j2^ , |2^ , p5] | . 
However, with the next-nearest-neighbor hopping term 
only a few results have been obtained with this technique 
p6f since the additional hopping amplitude exacerbates 
the sign problem (which exists even at half- filling) . 

In this paper, we present a systematic study of the evo- 
lution of the spectral function in the U — t — t' model at 
electronic densities ranging from half-filling ((n) = 1.0) 
to quarter- filling ((n) = 0.5) on a 6 x 6 lattice at an 
inverse temperature of (3t — 3(T = 0.33t). The pro- 
cedure we used to obtain these results is slightly differ- 
ent from those used before j2^j2(|, and we believe that 
it provides more details to the spectral function. First, 
approximately 100,000 QMC measurement sweeps were 
taken at each density (for a fixed T, t'/t, and U/t) to 
obtain accurate statistics for the imaginary time Green's 
functions. Once the spectral function was obtained from 
the standard ME procedure, we used this spectral func- 
tion as a seed for further ME analysis with a reduced 
coefficient of the entropy a plj , p7| . 

In order to work in the realistic strong coupling regime, 
sacrifices in the lattice size and temperature had to be 
made due to the sign problem |2S[ ]. The same set of 
measurements shown below in Sec. Ill could have been 
performed on an 8 x 8 lattice. However, it would 
have required a higher temperature ((3t = 2) to obtain 
good statistics on the imaginary time Green's functions 
thereby washing out parts of the quasiparticlc features 
that become prominent at lower temperatures. Since at 
half- filling and fit = 2, the 6x6 and 8x8 lattices give 
qualitatively similar results and we were able to reach 
a lower temperature using the smaller lattice, then we 
decided to use the 6x6 cluster throughout the paper. 
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FIG. 1. (a) Spectral functions ^(p,^) of the U — t — t' 
Hubbard model at U/t = 10, t'/t = -0.35, and T = rj/3, 
using QMC-ME techniques on a 6 x 6 cluster. The density 
is (n) = 1.0. From the bottom, the momenta are along the 
main diagonal from (0,0) to (ji,n), from there to (7r,0), and 
finally back to (0,0); (b) Energies of the dominant peaks in 
the spectral functions. The results in the PES region are ob- 
tained from a two-gaussian analysis of the QMC-ME results. 
The squares correspond to the incoherent part of the spec- 
trum (INC) and its error bars are not shown. The full circles 
form the quasiparticle (qp) band, and their diameter is pro- 
portional to the intensity. The error bars are given by the 
width at half the height of the gaussian corresponding to the 
q.p. peak at each momenta. The upper Hubbard band is also 
shown (open circles) without error bars; (c) density of states 
N(ui) obtained by summing the A(p,ui)'s. The inset shows 
the top of the q.p. peak in more detail. 
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III. QMC RESULTS 
A. Half-filling 

In Fig.l the QMC results obtained at half- filling for 
U = 10*, T = t/3, and t'/t = -0.35 are shown. Fig.la 
contains ^4(p, oj) after the ME analysis of the QMC data. 
The asymmetry of most of the dominant peaks suggests 
that they are a combination of at least two features lo- 
cated at different energies. This is reasonable since pre- 
vious calculations for the U — t model (see e.g. Ref. 
p5[) have established that below the chemical poten- 
tial A(p,u) is made out of a quasiparticle (q.p.) peak 
near /x, carrying a small fraction of the total weight, and 
a broad incoherent (INC) feature at larger binding en- 
ergies containing the rest of the weight. A(p,u>) have 
been fitted by two gaussians with positions, widths and 
weights adjusted to match the ME result, keeping the 
overall weight constant. The same study was performed 
above the chemical potential in the inverse photoemission 
(IPES) regime but such analysis did not provide interest- 
ing information since the weight of the peak the closest 
to n is very small, specially away from half-filling. Then, 
within the accuracy of the present QMC-ME study, the 
upper Hubbard band (UHB) is described as just contain- 
ing a broad featureless peak for each momentum. The 
two-peak decomposition analysis reported here will, thus, 
be limited to the PES part of the spectrum. However, 
note that at least at half-filling it is likely that the IPES 
portion of the spectrum contains a (low intensity) q.p.- 
like feature. In particular as t'/t — > particle- hole sym- 
metry must be recovered, and here the peaks appearing 
in the PES region must also exist in the IPES regime. 
Nevertheless, in Sec. III.B below we will show that at 
finite hole densities the q.p. band at the bottom of the 
UHB carries such a small weight that it can be neglected. 

Following this fitting procedure, in Fig. lb the energy 
position of the peaks as a function of momentum is 
shown. As anticipated, below \i a feature at the top 
of the valence band appears. It is natural to associate 
this peak with a "quasiparticle band" . The bandwidth is 
roughly 2t, i.e. much smaller than for free electrons on a 
lattice. However, it is not as small as predicted by t — J 
model calculations jf| , which may be due to the influence 
of t' or a finite U/t. The top of the q.p. band is located 
along the main diagonal from (0,0) to (tt, tt), in agree- 
ment with previous literature ||. The influence of the 
nonzero next-nearest-neighbor (NNN) hopping t' appears 
in the energy position of the q.p. band at (tt, 0) which is 
deeper in energy than (27r/3, 2ir/3) and (n/3, tt/3) (while 
for t'/t — 0.0, (tt, 0) is located very close to the top of 
the band in contradiction with ARPES experiments for 
cuprate insulators Q ) . It is likely that on larger clusters 
(tt/2, tt/2) would correspond to the actual top of the q.p. 
band, as predicted before [[|. The overall shape of the 
q.p. band is similar to that calculated in the t — J model 
using the self-consistent Born approximation |l(J . 
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FIG. 2. (a) Spectral functions A(p, to) of the 17 — t — t' 
Hubbard model at U/t = 10, t'/t = -0.35, and T = t/3, 
using QMC-ME techniques on a 6 x 6 cluster. The density 
is (n) = 0.94. From the bottom, the momenta are along the 
main diagonal from (0,0) to (tt, tt), from there to (tt,0), and 
finally back to (0,0); (b) Energies of the dominant peaks in 
the spectral functions. The results in the PES region are ob- 
tained from a two-gaussian analysis of the QMC-ME results. 
The squares correspond to the incoherent part of the spec- 
trum (INC) and its error bars are not shown. The full circles 
form the quasiparticle (qp) band, and their diameter is pro- 
portional to the intensity. The error bars are given by the 
width at half the height of the gaussian corresponding to the 
q.p. peak at each momenta. The upper Hubbard band is also 
shown (open circles) without error bars; (c) density of states 
N(u>) obtained by summing the A(p,uj)'s. 
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It is interesting to note that the q.p. band includes 
momentum (n,n), although with a very small intensity. 
This is reasonable for two reasons: (i) first, the presence 
of a charge gap in the spectrum, as clearly seen in Fig. lb, 
suggests that bands starting, say, at (0, 0) cannot simply 
disappear at some other momentum by crossing fi. By 
continuity they have to extend all along the Brillouin 
zone (unless their q.p. weight vanishes); (ii) in addition, 
the presence of strong antiferromagnetic correlations im- 
plies a doubling of the unit cell of dynamical origin and 
thus "shadow band" features should appear in the spec- 
trum, as discussed in several models with strong anti- 
ferromagnetic correlations (l| jl^,|l|,|l|,||,|oU29|j3ltl . In 
other words, AF introduces an extra symmetry in the 
problem that links p with p + (tt, tt). These shadow fea- 
tures are certainly weak compared to the rest of the band, 
but are nevertheless present in the spectrum of Fig. lb. 

The rest of the spectral weight is contained in the 
broad incoherent feature at energies deeper than the q.p. 
peak, and also in the UHB. The latter presents some 
nontrivial momentum dependence near (tt, tt). This be- 
havior seems a remnant of the results for noninteract- 
ing electrons. Finally, Fig.lc contains the density of 
states (DOS) N(uj) = £ p A(p,w) at half-filling. The 
gap A ~ 5t and the sharp q.p. band can be clearly iden- 
tified. 
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B. Finite Hole Density 

Fig. 2 contains the QMC-ME results obtained at den- 
sity (n) = 0.94, with the rest of the parameters as in 
Fig.l. In Fig. 2a the spectral weight A(p,u>) is shown. 
There are several interesting features in this result. For 
example, the upper Hubbard band has lost weight com- 
pared with results at (n) = 1.0. This is not surprising 
since as the hole density grows the chances that an elec- 
tron injected in the lattice will populate already occupied 
sites is reduced. 

The most interesting physical consequences of hole 
doping are obtained in the PES region. Here the chemi- 
cal potential has moved to the top of the valence band, 
more specifically into the q.p. band. Fig. 2b shows the re- 
sults of a two-peak analysis of the PES region similar to 
that performed at half-filling. Interpolating results along 
the main diagonal it is observed that the q.p. peak at 
momentum p = (7r/2,7r/2) should be approximately at 
the chemical potential for this density, while (tt, 0) is still 
about 2t below. The q.p. band has distorted its shape 
in such a way that considerable weight has been moved 
above /i in the vicinity of (tt, tt). Actually, now the q.p. 
band resembles a more standard tight-binding dispersion 
although with a renormalized hopping amplitude smaller 
than the bare one. Note, however, that the position of 
the q.p. peak at momentum (0, 0) is somewhat patholog- 
ical since it is located above the chemical potential. We 
do not have an explanation for this anomalous behavior, 
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FIG. 3. (a) Spectral functions ^(p,^) of the U — t — t' 
Hubbard model at U/t = 10, t'/t = -0.35, and T = t/3, 
using QMC-ME techniques on a 6 x 6 cluster. The density 
is (n) = 0.84. From the bottom, the momenta are along the 
main diagonal from (0,0) to (tt,tv), from there to (tt,0), and 
finally back to (0,0); (b) Energies of the dominant peaks in 
the spectral functions. The results in the PES region are ob- 
tained from a two-gaussian analysis of the QMC-ME results. 
The squares correspond to the incoherent part of the spec- 
trum (INC) and its error bars are not shown. The full circles 
form the quasiparticle (qp) band, and their diameter is pro- 
portional to the intensity. The error bars are given by the 
width at half the height of the gaussian corresponding to the 
q.p. peak at each momenta. The upper Hubbard band is also 
shown (open circles) without error bars; (c) density of states 
N(ui) obtained by summing the A(p,uj)'s. 
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FIG. 4. (a) Spectral functions A(p, w) of the U — t — t' 
Hubbard model at LT/t = 10, /j'/t = -0.35, and T = i/3, 
using QMC-ME techniques on a 6 x 6 cluster. The density 
is (n) = 0.78. From the bottom, the momenta are along the 
main diagonal from (0,0) to (ji,n), from there to (7r,0), and 
finally back to (0,0); (b) Energies of the dominant peaks in 
the spectral functions. The results in the PES region are ob- 
tained from a two-gaussian analysis of the QMC-ME results. 
The squares correspond to the incoherent part of the spec- 
trum (INC) and its error bars are not shown. The full circles 
form the quasiparticle (qp) band, and their diameter is pro- 
portional to the intensity. The error bars are given by the 
width at half the height of the gaussian corresponding to the 
q.p. peak at each momenta. The upper Hubbard band is also 
shown (open circles) without error bars; (c) density of states 
N(uS) obtained by summing the A(p,w)'s. 
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FIG. 5. Density of states N(u>) corresponding to the 
U - t - t' Hubbard model at U/t = 10, t'/t = -0.35, and 
T — t/3, using QMC-ME techniques on a 6 x 6 cluster. The 
density is (n) = 0.65. 

which likely is caused by the Maximum Entropy proce- 
dure. The incoherent part of the spectrum (INC) is also 
clearly visible in the calculation. At this density it re- 
mains entirely filled, i.e. the electrons removed from the 
system have been taken from the q.p. band. Finally, note 
that the shadow features are no longer prominent. This 
is correlated with a substantial reduction of the antiferro- 
magnetic correlation length £af at this density and tem- 
perature compared with the results at half-filling. This is 
not in contradiction with claims that recently observed 
features in ARPES are induced by antiferromagnetism 
|]r7| , ^8| since at very low temperatures £af likely remains 
robust for a larger hole density region near half-filling 
than it occurs at the relatively high temperature T = t/2>. 
In Fig.2c the DOS for (n) = 0.94 is shown. The q.p. band 
is sharp, and \i lies slightly to the right of the top of this 
band (see inset of Fig. 2c). A pseudogap generated by the 
large coupling U/t is still clearly visible. 

Fig. 3a contains the QMC-ME results obtained at (n) = 
0.84. Compared with the results at (n) = 0.94, now /j, 
is located deeper inside the q.p. band according to the 
two-gaussians analysis presented in Fig. 3b. Interpola- 
tions clearly show that p = (tt/2,tt/2) for this band is 
now in the IPES region, while (tt, 0) is at or very close 
to fi but still in the PES regime. The incoherent weight 
in PES remains robust. Fig. 3c shows the DOS at this 
density, fi is located slightly to the left of the q.p. band 
maximum, the pseudogap has virtually melted, and there 
is only a tiny trace of the upper Hubbard band. The 
presence of (tt, 0) very close to /i is correlated with ex- 
perimental results showing "flat bands" near the Fermi 
energy at optimal doping R p^ , ^flj33| . 

Fig. 4 contains the QMC results obtained at (n) = 0.78. 
This density is particularly interesting since now seems 
to have crossed most of the q.p. band [Q. This can be 
observed directly in the spectral weight (Fig. 4a), and in 
its two-peak analyzed output (Fig. 4b). This detail will 
have interesting consequences for some theories of high- 
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FIG. 6. (a) Spectral functions A(p, uj) of the U — t — t 
Hubbard model at U/t = 10, t'/t = -0.35, and T = t/3, 
using QMC-ME techniques on a 6 x 6 cluster. The density 
is (n) = 0.47. From the bottom, the momenta are along the 
main diagonal from (0,0) to (tt,tv), from there to (7r,0), and 
finally back to (0,0); (b) Energies of the dominant peaks in 
the spectral functions. The solid line is the dispersion in the 
noninteracting limit U/t — 0.0; (c) density of states N(u) 
obtained by summing the A(p,u)'s. 



T c uprates, as described later in this paper. The DOS 
is shown in Fig. 4c. The same trend is also observed 
at (n) — 0.65 for which here only the DOS is shown 
(Fig. 5). At this density the q.p. band is empty and the 
INC becomes sharper presumably due to its proximity to 
the chemical potential. It is interesting to note that the 
chemical potential in Fig. 5 is located at a sharp minimum 
in the density of states. 

Results at (n) = 0.47 (Fig. 6) show that the q.p. band 
has melted and at this low density the spectral weight 
is dominated by features that can be traced back to 
the incoherent part of the PES spectrum at half-filling 
(Fig. 6a). In Fig. 6b the position of the dominant peak 
is given as a function of momentum, and in Fig. 6c the 
DOS is shown. It is interesting that now the results are 
nicely fitted by the non-interacting tight-binding disper- 
sion (solid line in Fig. 6b). Thus, the "free electrons" limit 
is approximately recovered in the strong coupling regime 
at "quarter-filling" . Again, note that the peak structure 
in the spectrum seems to have emerged from the inco- 
herent part of the valence band observed at half-filling. 



IV. COMPARISON OF QMC-ME RESULTS WITH 
ARPES EXPERIMENTS 

As discussed in the Introduction, ARPES results have 
shown that the quasiparticle peak at p = (n, 0) is very 
sensitive to hole doping at least in the underdoped regime 
For the same densities, the q.p. dispersion along 
the main diagonal from (0, 0) to (it, tt) does not change 
as much. In this section it will be investigated if the nu- 
merical results of Sec. Ill are compatible with the ARPES 
data. 

In Fig. 7, we schematically show what occurs upon dop- 
ing when a tight binding dispersion including a NN hop- 
ping amplitude t* , and a NNN amplitude t'* with a ratio 
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FIG. 7. Density evolution of a non-interacting band, re- 
ferred to the chemical potential, as its density changes (shown 
in the figure). For details see the text. 
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FIG. 8. (a) Energy of the quasiparticle (referred to the 
chemical potential) in units of t for the U — t — t' Hubbard 
model using the parameters of Figs. 1-6. The results here cor- 
respond to momentum p = (tv,0). The solid circles, squares 
and open circles are the q.p., INC and UBH positions for this 
momentum; (b) quasiparticle spectral function corresponding 
to momentum p = (tt, 0) for several densities (n) (indicated). 
The small triangles show the positions of the dominant peaks 
in the two-gaussian analysis of the PES data. 
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FIG. 9. Energy of the quasiparticle (referred to the chemi- 
cal potential) in units of t for the U — t — t 1 Hubbard model 
using the parameters of Figs. 1-6. The results here correspond 
to momentum p = (2tt/3, 2tt/3). The solid circles, squares 
and open circles are the q.p., INC and UBH positions for this 
momentum; (b) quasiparticle spectral function corresponding 
to momentum p = (27r/3, 27r/3) for several densities (n) (in- 
dicated). The small triangles show the positions of the dom- 
inant peaks in the two-gaussian analysis of the PES data. 

t'*/t* = —0.35 is used. The reason for this exercise is 
that the QMC-ME results of Sec. Ill have shown that 
rapidly upon hole doping the quasiparticle dispersion re- 
sembles that of a renormalized non- interacting set of elec- 
trons. Then, it is interesting to analyze what would occur 
with such dispersion as the density within the q.p. band 
changes (i.e. the density (n) shown in Fig. 7 must be con- 
sidered as the filling of the q.p. band, rather than the 
density of the whole system). It is clear from the fig- 
ure that (7r, 0) is intrinsically more sensitive to doping 
than the main diagonal. Actually the results resemble in 
part those of Marshall et al. specially the experimental 
ARPES dispersion corresponding to slightly underdoped 
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FIG. 10. Best fits of the QMC-ME quasiparticle dispersion 
corresponding to densities (n) — 0.84 and (n) = 0.78. The 
results at (n) = 1.0 are taken from Ref. [10] arbitrarily locating 
the top of the band at the chemical potential. 

samples of Bi2212 with T c = 85K, as well as previ- 
ous experiments for optimally doped samples 0, where 
the q.p. band crosses the chemical potential between 
(tt, tt) and (w, 0). This lead us to conjecture that the 
ARPES experiments in the slightly underdoped regime 
may have observed the hole filling of a "free- like", but 
narrow, q.p. band. If this speculation is correct, then 
it is here predicted that (tt, 0) will eventually cross the 
chemical potential when the system becomes overdoped, 
as Fig. 7 shows. To the best of our knowledge there are 
no ARPES results in this regime. Experimental work for 
overdoped cuprates would clarify the issue of whether the 
"flat bands" found in the optimal regime remain locked 
near \i or smoothly cross the Fermi energy as the hole 
doping increases. Our results favor the latter, at least 
within the resolution of the QMC-ME methods. 

Further evidence that (tt, 0) travels across the chem- 
ical potential is given in Fig. 8a, b. There A(p,u>) with 
p = (tt, 0), taken from Figs.la-6a, is shown again adding 
also results for densities (n) = 0.89 and 0.65 for com- 
pleteness. The peak at this momentum belonging to the 
q.p. band crosses fi at a density (n) « 0.82. Then, 
as the hole density grows a movement up in energy of 
the p = (tt, 0) quasiparticle is clearly observed, as in 
ARPES experiments j3f|. At a larger hole density, such 
as quarter-filling, the q.p. peak, now in the IPES regime, 
substantially reduces its weight and disappears. At this 
density the feature associated with the incoherent part 
of the spectrum at half-filling is now located close to the 
Fermi energy, contributing to the dispersion that resem- 
bles the U = result (Fig. 6b). The UHB rapidly looses 
weight moving away from half-filling. Fig. 9a, b show sim- 
ilar results but for p = (27r/3, 27r/3), which is represen- 
tative of the behavior along the main diagonal on the 
6x6 cluster. The q.p. associated to this momentum 
crosses the chemical potential at a smaller hole density 
((n) w 0.95). 

Summarizing the results of this section, in Fig. 10 three 
representative q.p. dispersions obtained with QMC-ME 
methods (Sec. Ill) are shown. For the half-filled case 
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the results are accurately known from the self-consistent 
Born approximation [[nj, and (tt,0) is clearly below the 
chemical potential, which here is arbitrarily located at 
the top of the band. Note that the intensity of the q.p. 
peak is much weaker near (tt, tt) than for other momenta 
at half-filling. For small hole doping densities, such as 
(n) = 0.84, a "free-like" but narrow quasiparticle band 
is observed in QMC simulations (also in studies of the 
U — t model J36| , p5[ ) and (tt, 0) is now at the Fermi en- 
ergy. The vicinity of (tt, 7t) is the region affected the 
most by hole doping due to the reduction of antiferro- 
magnetic correlations. Finally, at lower densities, such 
as (n) = 0.78, the chemical potential has crossed the 
q.p. band and now (tt, 0) is above the Fermi energy. We 
consider that the scenario depicted in Fig. 10 for the den- 
sity evolution of the q.p. dispersion can provide a sim- 
ple qualitative explanation for the ARPES results found 
experimentally. It predicts that the q.p. peak at (tt, 0) 
will eventually disappear from the ARPES signal into the 
IPES region, an effect that could be tested experimen- 
tally. At densities intermediate between (n) = 1.0 and 
0.84, the peak at (tt, 0) should evolve smoothly between 
these two limits, and thus close to half-filling "pseudo- 
gap" features will likely appear in simulations if lower 
temperatures could be reached. Reducing the tempera- 
ture below T = t/3 is particularly important for the gen- 
eration of strong spin fluctuations which are crucial to 
induce the quasiparticle dispersion at half-filling shown 
in Fig. 10. Note that "shadow" features should actually 
be present at (n) = 0.84 but they are too weak to be 
detected by QMC-ME methods, and probably also by 
ARPES experiments which have large backgrounds in 
their signals. Thus, the existence of hole pockets can- 
not be shown from the current Monte Carlo simulations 
available. 



V. IMPLICATIONS OF QMC RESULTS FOR 
SOME THEORIES OF HIGH-T C 

In this section the implications of our QMC results for 
recently proposed scenarios for high-T c cuprates will be 
discussed. Two features found in the simulation will be 
important, namely the "tight-binding" shape of the q.p. 
band at finite hole density (although with renormalized 
parameters) and the crossing of the q.p. band by fi for 
densities in the approximate range 0.70 < (n) < 1.0. 



A. Large peak in the DOS 

As explained in the introduction, the presence of "flat" 
regions in the experimental normal-state q.p. dispersion 
is a remarkable feature of the phenomenology of hole- 
doped cuprates W. These flat bands are located around 
momenta p = (tt, 0) and (0, tt), and at optimal doping 
they are ~ lOmeV below the Fermi energy M. Studies 



of holes in the 2D t — J and Hubbard models at and 
away from half-filling have suggested that antiferromag- 
netic correlations may play an important role in the gen- 
eration of these features |3^j3fj| , p3| , p^6| . Our results for 
the U — t — t' model suggest the presence of flat bands 
in the regime close to (n) = 0.84, although care must 
be taken with the finite resolution of the ME-generated 
q.p. peaks. But even if these flat regions were not quan- 
titatively described by one band electronic models, the 
intrinsic small bandwidth of the q.p. band at small hole 
density could be enough to induce a large peak in the 
DOS which can be used to enhance the superconducting 
critical temperature, once a source of hole attraction is 
identified. This leads to a natural explanation for the ex- 
istence of an "optimal doping" which in this framework 
occurs when the peak in the q.p. band is reached by 
/i It is also natural to label as "underdoped" the 

regime where /i is to the right of the "flat band" peak 
(i.e. at higher energies), and "overdoped" when it is to 
the left (i.e. at lower energies) 
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FIG. 11. (a) N(lu) for the 2D t - J model obtained with 
exact diagonalization techniques averaging results for clusters 
with 16 and 18 sites, at J ft — 0.4 and for the densities indi- 
cated. The (^-functions were given a width r\ — 0.25t. Similar 
results were found at other values of J/t; (b) N(ui) for the 
one band Hubbard model obtained on a 4 x 4 cluster using 
QMC and ME, but without reducing the coefficient of the en- 
tropy as in Sec. III. The temperature is T = i/4 and U/t = 12. 
Densities are indicated. 

Previous analytical calculations jl2| assumed the sur- 
vival of the large DOS obtained at half-filling as the 
density of holes grows. The results of our simulations 
(Sec. Ill) allow us to study the evolution of the DOS with 
hole density and judge if these ideas are realistic. The 
presence of a robust q.p. peak is certainly confirmed by 
our results, and Figs.lb-4b show that as (n) decreases 
from 1, the q.p. band is crossed by /i as conjectured be- 
fore @. The DOS of the standard t — J and Hubbard 
models (with t' = 0) also have a large peak in the DOS, as 
found in previous numerical simulations. For complete- 







ness, here those results are also presented. In Fig. 11a, 
N(u>) for the 2D t — J model obtained with exact diago- 
nalization (ED) techniques is shown at several densities 
||. At half-filling, a large DOS peak appears at the top 
of the valence band [B2|. Note that substantial weight 
exists at energies far from u, i.e. the large peak carries 
only a fraction of the total weight, in agreement with 
the QMC simulations of Sec. III. In the t — J model, the 
maximum in the DOS is not strictly located at the top 
of the valence band but at slightly smaller energies j52| . 
This effect is enhanced by adding t' hopping terms to the 
t — J Hamiltonian, as shown in the simulation results of 
Sec. III. As (n) decreases, the peak in Fig. 11a is now much 
broader but it remains well-defined. At (n) ~ 0.88, u is 
located close to the energy where N(w) is maximized. At 
(n) ~ 0.75, fi moves to the left of the peak. These results 
are quantitatively similar to those of our QMC simula- 
tions for the U — t — t' model, and also with results for 
the t'/t = 0.0 Hubbard model (Fig. lib). Thus, while ac- 
curate extrapolations to the bulk limit are difficult, the 
simple qualitative picture emerging from these studies, 
namely that strong correlations generate a q.p. peak in 
the DOS which is crossed by fj,, seems robust. 

As explained before, this q.p. -peak-crossing of the DOS 
is important since once a source of hole attraction exists 
in the system, a superconducting (SC) gap would open 
at /i, and the resulting T c could be enhanced due to the 
large number of states available. The numerical results, 
both QMC and ED, are thus compatible with scenar- 
ios where a large T c is obtained due to an increase in 
N(p). Since the peak width increases substantially with 
hole density, strictly speaking the rigid band filling of 
the half-filled hole dispersion is invalid. However, such 
an approximation seemed to have captured part of the 
qualitative physics of the problem, since the DOS peak 
is not washed out by a finite hole density. 

In the cases discussed here, i.e. the U — t — t' and t — J 
models, calculations of the spin-spin correlations show 
that £af is approximately a couple of lattice spacings 
when fi is located near the DOS peak, becoming smaller 
as the overdoped regime is reached. Thus, a nonzero £af 
and n near a large DOS peak are correlated features. It 
is in this respect that scenarios where AF correlations 
produce a large DOS that enhances T c |l2| differ from 
vH theories where divergences in the DOS are caused 
by band effects already present before interactions are 
switch on [l^j . 

While the existence of a robust peak in the DOS is in 
good agreement with ARPES data it is in apparent 
disagreement with specific heat studies for YBCO [p8| . 
The lack of p-resolution in the specific heat measure- 
ments may solve this puzzle. Actually, angle-integrated 
PES results for Bi2212 do not show the sharp flat features 
found in ARPES for the same material |p9| , p0[ . Similar 
effects may affect the specific heat data which should be 
reanalyzed to search for DOS large peaks. 



B. Kondo resonances vs AF induced quasiparticle 

Previous QMC-ME studies of the Hubbard model for 
t'/t = 0.0 [[36 41 T3] reported results qualitatively similar 
to those shown in Fig. lib, where the DOS obtained on 
a 4 x 4 cluster is presented for the one band Hubbard 
model at U/t = 12, and T = i/4. The ME technique used 
in Fig. lib is the same as in those previous simulations, 
and it does not have the resolution of the present ME 
variation discussed in Sec. III. The crossing of a peak in 
the DOS by /i is clearly observed in these ME simulations 
but at half- filling there are no q.p. peaks and Fig. lib can 
thus be naively interpreted as the "generation" by doping 
of a Kondo- like peak at the top of the valence band which 
does not exist at (n) = 1. However, our current results 
(Fig.lc) actually show that at (n) = 1.0 a well-defined 
q.p. peak is present, as predicted by a variety of studies 
of the zero temperature t — J model with one hole ||^5) ■ 
Also experimentally in the cuprates it has been already 
established j3{| that the states observed in PES upon 
doping are already present in the insulator and are not 
Kondo resonances. Thus, Figs.lb-4b provide evidence 
that the q.p. features observed at (n) < 1 are smoothly 
connected to peaks already present at half-filling. Similar 
conclusions have also been obtained in geometries other 
than the 2D square lattice, such as a t — J ladder [BjlEgl, 



C. d x 2_ y 2 in the Hubbard model 

The results of Figs.lb-6b show that the q.p. band is 
very sensitive to doping, at least at the temperatures and 
couplings used in this QMC-ME study. In particular a 
few percent hole doping is enough to transform the half- 
filling dispersion, containing the extra symmetry induced 
by long-range AF order, into a tight-binding-like disper- 
sion although with a small hopping amplitude. These 
results are similar to others previously published in the 
literature. For example in Fig. 12a, b the q.p. dispersion 
at (n) = 1, and at (n) ~ 0.87 for the t — J model is repro- 
duced from Ref. |25| (for Hubbard model results, see Ref. 
p6| , ^3[ ). Upon doping, vestiges of the flat regions remain, 
inducing a large peak in the DOS (see Fig. 11). However, 
the region around p = (-7T, it) has changed substantially, 
i.e. the AF shadow region observed at (n) = 1 reduced 
its intensity and considerable weight was transferred to 
the IPES region. The q.p. dispersion at (n) ~ 0.87 of the 
t — J model can also be fitted by a tight-binding nearest- 
neighbors (NN) dispersion with a small effective hopping 
likely associated with J. Then, both previous literature 
results and the present simulations agree on the qualita- 
tive aspects of the evolution with doping of A(p,w). 

The changes in the q.p. dispersion with hole doping 
can be interpreted in two ways: First, note that A(p, lu) is 
influenced by matrix elements of bare fermionic operators 
connecting states with N and N ± 1 particles. This is 
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FIG. 12. (a) q.p. energy vs momentum obtained at 
half-filling using the t — t' — J (from Refs. [12,32]). The re- 
sult shown, that we call 6af(p), is a good fit of Monte Carlo 
data on a 12 x 12 cluster at J/t = 0.4; (b) q.p. dispersion 
vs momentum at (n) = 0.87 and J/t = 0.4 using exact di- 
agonalization of 16 and 18 sites clusters for the t — J model 
(from Ref.[25]). The open (full) circles are IPES (PES) re- 
sults. Their size is proportional to the peak intensity. The 
solid line is the fit ejvjv(p) described in the text; (c) V(x) 
along the x-axis after Fourier transforming the smeared po- 
tential V(p) = <5(p — Q) (see text). £af is given in lattice 
units. 

important when the q.p. weight is small, i.e. when the 
state c pa -\gs)N does not have a large overlap with the 
ground state of the N — 1 particles subspace c/s) a^— l 
(\gs)N being the ground state with TV particles). If the 
hole excitation is instead created by a new operator 7 pcr 
that incorporates the dressing of the hole by spin fluctu- 
ations, then 7 P o-|gs) n may now have a large overlap with 
|<?s)iv-i @,(f5|- In other words, if the dressed hole state 
resembles an extended spin polaron, then the physics de- 
duced from PES studies, which rely on the sudden re- 
moval of a bare electron from the system, may be mis- 
leading. To the extent that spin polarons remain well- 
defined at finite density, the use of 7 p(T will induce spec- 
tral weight rearrangements between the PES and IPES 
regions, and the results at small hole density could re- 
semble those at half-filling after such spectral weight re- 
distribution takes place. If this idea were correct, then 
the study of superconductivity and transport, both regu- 
lated by dressed quasiparticles, could indeed be handled 
by filling a rigid band given by caf(p)- K is likely that 
this idea will work in the underdoped regime where £af 
is robust, and results for 2D clusters |Q and ladders 
p6[ already support these claims. Note that the fact 
that the bandwidth at finite hole density remains much 
smaller than 8t shows that strong correlations still play 
an important role in the dispersion. 

However, an alternative is that the difference between 
Figs. 12a and 12b for the t — J model and Figs. lb and 2b 
for the U — t — t' model correspond to an intrinsic change 
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FIG. 13. (a) T c vs [i obtained with the BCS gap equation, 
with [i = as the chemical potential corresponding to the 
saddle-point. The dashed line corresponds to results using 
caf(p), and the solid line to ejvjv(p). Interpolations between 
these two extreme cases can be easily constructed using two 
band models with weights regulated by Z-factors (actually 
we calculated T c using a two-pole Green's function in the gap 
equation and the result smoothly interpolated between those 
shown in the figure). Note the presence of both d x 2 _ y 2 -wave 
and extended s-wave SC; (b) T c for d-wave SC vs the percent- 
age x of filling of the q.p. band (i.e. not of the full system). 
The dashed line is the same as in (a). The solid line corre- 
sponds to the BCS gap equation result making zero the weight 
Z p of states in the q.p. dispersion that are at energies from 
the saddle point larger than 2.5% of the total bandwidth (i.e. 
basically including states only in a window of energy ~ 125K 
around the flat regions). The dot-dashed line is the same but 
using a window of ~ 250 AT around the flat regions. 

in the q.p. dispersion as (;af decreases. This is the less 
favorable case for the real-space pairing approaches 
which are constructed at half-filling, and thus we should 
analyze this possibility in detail here. For this purpose, 
we have applied the standard BCS formalism to an ef- 
fective model with a low density of quasiparticles hav- 
ing a dispersion ejvAf(p)/eI / = — 0.2(cosp x + cosp y ) (as- 
suming t = OAeV), which roughly reproduces the dom- 
inant features of Fig. 12b (the reader should not confuse 
this dispersion with that of free particles whose band- 
width is four times larger. A large amount of incoher- 
ent weight is still observable at the density of Fig. 12b 
although it is not shown explicitly. The same occurs 
in Fig. 2b for the U — t — t' model). As described be- 
fore, to study T c we should include a NN attraction in- 
duced by AF between these quasiparticles. While naively 
it may seem dubious to use the same interaction both 
at and away from half-filling, the hole-hole potential 
should not be much affected at distances shorter than 
£af- This can be illustrated by a real space analysis 
(Fig. 12c) of a smeared (5-function potential of AF origin 
V(q) — £af/[1 + £,af( < 1~ Q) 2 ]j where the lattice spacing 
is set to 1, and Q = (7r,7r). Fig. 12c shows that the NN 
potential (x = 1) does not change noticeably as £af is 
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reduced, while V(x > 1) is rapidly suppressed. Then, us- 
ing the same NN form of the potential for many densities 
should not be a bad approximation J47] . 

To analyze the stability of previous calculations let us 
then use the tight-binding dispersion £nn(p) obtained 
from the numerical analysis in the BCS gap equation. 
Solving numerically the gap equation, T c is shown in 
Fig. 13a. As reported before |12[, superconductivity at 
T c ~ 80 — lOOivf in the d x 2_ y 2 channel appears natu- 
rally if the hole dispersion eaf(p) at half- filling is used. 
Nevertheless, if instead €nn(p) is used, the vestiges of 
the flat bands present in this narrow dispersion produce 
T c ~ 30K which is still large pq |. Even more remark- 
able is the fact that the d x 2_ y 2 character of the SC state 
is maintained. This result can be understood noticing 
that a combination of e nn (p) with an attractive NN po- 
tential effectively locates the Hamiltonian in the family 
of "t — U — V" models with U repulsive and V attrac- 
tive, where it is known that for a "half-filled" electronic 
band the dominant SC state is d x 2_ y 2-wave p9| . In other 
words, when \i is at the flat region of the ejvjv(p) disper- 
sion it approximately corresponds to an effective "half- 
filled" q.p. band resembling a free electron dispersion 
(but with smaller bandwidth), leading to a d x 2_ y 2-w&ve 
SC state (Fig. 13a). Then, even if the q.p. dispersion 
changes substantially with doping near the Q = (it, tt) 
point, such an effect does not seem to alter the main 
qualitative features found in previous studies [fl2|| . This 
is a remarkable result and it is caused by the rapid de- 
population of the q.p. band as the hole density grows. 
Then, while (n) governs the actual global density of the 
system, there is a subtle hidden density (i.e. the popu- 
lation of the q.p. band) that may influence considerably 
on the physics of the problem. At (n) = 0.78 and for the 
U—t—t' model, this q.p. band is nearly empty while glob- 
ally the system is still close to half-filling. While the use 
of £af(p) or £nn{p) is certainly a rough approximation, 
the existence of a robust T c and d-wave superconductiv- 
ity in both cases suggests that similar results would be 
obtained using more realistic models for the hole pairing 
interaction. 



D. Influence of shadow bands on SC 

The present analysis also shows that the AF "shadow" 
regions of the half-filling hole dispersion caf(p) are not 
crucial for the success of the real-space approach. Us- 
ing cnn(p), which does not contain weight in PES near 
(ir, 7r), T c is still robust and the d-wave state remains 
stable, as shown in the previous subsection. To establish 
this result more clearly, we analyzed T c using €af(p) but 
modulating the contribution of each momentum with a p- 
dependent weight Z p in the one particle Green's function. 
We considered the special case where Z p is zero away 
from a window of total width W centered at the saddle 
point, which is located in the flat bands region. Inside 



the window W, the weight is maximum i.e. Z p = 1. Such 
a calculation also addresses indirectly possible concerns 
associated with widths ~ (eaf(p) — £_f) 2 that q.p. peaks 
would acquire away from half-filling in standard Fermi 
liquids (ep is the Fermi energy). Results are shown in 
Fig. 13b, for d-wave SC. Note that even in the case where 
W is as small as just 5% of the total bandwidth (itself 
already small of order 2J), T c remains robust and close 
to 70K. Then, it is clear that the dominant contribu- 
tion to T c comes from the flat regions and the shape 
of the q.p. dispersion away from them has a secondary 
importance for the success of the real-space approach. 
However, it is important to remark that the calculation 
described in this section showing that shadow features 
are not much important for the actual value of T c does 
not mean that AF correlations can be neglected. On the 
contrary, the whole Antiferromagnetic Van Hove scenario 
Jl2| and other similar approaches are based on the notion 
that pairing is caused by spin fluctuations |30) . 

E. SC in the overdoped regime 

Finally, novel predictions obtained in the regime where 
the q.p. band is almost fully crossed by [i (as observed in 
Sec. Ill for (n) between 0.78 and 0.66) are here discussed. 
From the point of view of the q.p. band this regime is 
"dilute" but, again, this should not be confused with the 
bottom of the whole spectrum since a large amount of 
weight remains in the incoherent part of A(p, to). At this 
density a standard BCS gap equation analysis applied to 
a q.p. dispersion either constructed at half-filling, as in 
Ref. p3 , or phenomenologically obtained from our data 
(ejvjv(p)), and supplemented by nearest-neighbor attrac- 
tion induced by antiferromagnetism, shows that extended 
s-wave SC dominates over d x 2_ y 2-wave SC. This corre- 
sponds to the "overdoped" regime (Fig. 13a) i.e. to an 
overall density in the vicinity of (n) ~ 0.70 |pl| . This 
change in the symmetry of the SC state can be under- 
stood recalling once again that the tight-binding disper- 
sion of the q.p., even including renormalized amplitudes, 
supplemented by a NN attraction, formally corresponds 
to an effective H — U — V" model. It is well-known that 
in this model the SC state symmetry changes from d- to 
s-wave as the density is reduced away from half-filling to 
a nearly empty system p9| , p2[ . Actually the bound state 
of two particles on an otherwise empty lattice with a NN 
tight binding dispersion and NN attraction is s-wave. To 
the extend that the AF- or NN-dispersions survive up to 
~ 25% hole doping, as suggested by the numerical data 
of Sec. Ill as well as previous literature (Ref. |25| and ref- 
erences therein), scenarios based on the real-space inter- 
action of q.p.'s predict a competition between extended 
s-wave and d 2 .2_ J/ 2-wave SC in the overdoped regime. Re- 
cent calculations based on the analysis of the low elec- 
tronic density (n) <C 1 limit of the t — J model led to 
analogous conclusions [p3|. Our approach is based on a 
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very different formalism but it arrives to similar results, 
and, thus, a crossover from d- to s-wave dominated super- 
conductivity in overdoped cuprates could occur. Indeed 
recent ARPES data for overdoped Bi2212 have been in- 
terpreted as corresponding to a mixing of s- and d-wave 
states [[34}. More work should be devoted to this po- 
tential competition between d- and s-wave for overdoped 
cuprates. 

VI. SUMMARY 

In this paper the results of an extensive numerical 
study of the U — t — t' one band Hubbard model with 
t'/t = —0.35 have been presented. The regime of strong 
coupling U/t — 10 was analyzed. These parameters 
are fixed to reproduce ARPES data for the AF insu- 
lator. With maximum entropy techniques the spectral 
function A(p, lu) was studied for several electronic densi- 
ties. It was observed that as the hole density grows away 
from half-filling, the quasiparticle band acquires "non- 
interacting" features although with bandwidths substan- 
tially smaller than for the U = limit. As the hole 
density increases the q.p. peak at (tt, 0) rapidly changes 
its position relatively to the Fermi energy, in qualita- 
tive agreement with recent ARPES experiments (l],|2j . In 
the overdoped regime it is predicted that the flat bands 
should be crossed by the chemical potential and thus they 
should no longer be observed in ARPES studies. The 
present results also have implications for some theories 
of bigh-T c . The narrow q.p. band produces a large peak 
in the DOS induced by strong correlations, which sur- 
vives the presence of hole doping and it can be used to 
boost T c once a source of hole attraction is found. It 
is remarkable that both with the hole dispersion found 
at half- filling or the one observed using QMC-ME at, 
e.g., (n) — 0.84, superconductivity in the d x 2_ y 2 channel 
dominates. This result gives support to real-space pair- 
ing theories of high-T c |12|] , showing that the main ideas 
of the scenario are stable upon the introduction of a fi- 
nite hole density. The large DOS peak is crossed by the 
chemical potential as (n) is reduced. In the overdoped 
regime, i.e. when the q.p. band is almost empty, a possi- 
ble competition between extended-s and d x 2_ y 2 SC was 
discussed. 
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